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1 . 0 INTRODUCTION 


The  theoretical  analysis  of  exact  element  patterns  from  circum- 
ferential and  radial  slots  near  the  tip  of  a sharply  tipped  cone  has 
been  completed  and  presented  in  the  final  report  on  Contract  No.  N0019-73-C-0127 
(Bargeliotes , Kummer,  and  Villeneuve,  1974).  In  both  cases,  the  element 
patterns  have  been  obtained  by  summation  of  the  normal  mode  series 
provided  by  the  angular  eigenfunctions  of  the  boundary  value  problem 
expressed  in  spherical  coordinates.  Relatively  few  modes  are  required 
for  the  convergence  of  the  series  for  slot  locations  near  the  tip  of 
the  cone,  the  number  of  modes  increasing  with  increased 
distance  from  the  tip.  In  addition,  computational  difficulties  from 
accumulated  round-off  errors  also  hinder  the  applicability  of  the  modal 
series  program  for  element  locations  far  from  the  tip. 

The  normal  mode  series  may  be  derived  from  related  integrals 
evaluated  as  a sum  of  residues  taken  over  the  poles  of  the  integrands 
(Pridmo re-Brown,  1972).  This  integral  represer tation  of  the  fields 
lends  itself  to  asympotic  approximations  which  allow  the  separation 
of  the  diffraction  field  and  the  geometric  optical  field.  In  turn, 
the  diffracted  field  may  be  expressed  as  a product  of  some  frequency- 
independent  angular  functions  or  diffraction  coefficients,  om(8,  0o) , 
and  a function  f(r  + a)  / (ka)  Here  qd  is  the  exterior  cone 

angle,  a is  the  distance  from  the  tip  to  the  slot  measured  along  the 
generatrix,  and  k is  the  wave  number.  The  geometry  is  shown  in  Figure 
1.  For  slots  far  from  the  tip  this  representation  of  the  diffraction 
field  is  better  suited  for  numerical  computation  than  the  modal  series. 
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a 


The  approximate  evaluation  of  the  contour  integrals  representing 


the  0-component  of  the  electric  field  of  a circumferential  slot 
resulted  in  explicit  expressions  of  the  field  in  terms  of  the  diffraction 
coefficients  mentioned  above  (Pridmore-Brown , 1972).  An  attempt  to 
utilize  these  expressions  for  the  numerical  computation  of  the  diffraction 
coefficients  and  the  field  itself  revealed  some  inconsistancies  or 
omissions  among  the  given  parameters.  In  an  effort  to  correct  these 
deficiencies,  the  complete  derivation  of  the  expressions  was  undertaken. 
Additionally,  for  a complete  pattern  analysis  or  synthesis  technique, 
similar  expressions  are  also  required  for  the  diffraction  coefficient 
associated  with  the  ({(-polarization . 

The  integral  expressions  for  the  scalar  potential  functions  that 
determine  the  radiation  fields  due  to  a cosm<{>  excitation  are  of  the 
form  (Pridmore-Brown,  1972) 


/- 


vdv 


A"1  (0,0O)  i-%  (ka) 


2 1 
V - 4 


where  the  angular  factor  A ™ (6,&0)  is  a ratio  of  associated  Legendre 
functions  and  the  integrals  are  taken  over  a hairpin  contour  enclosing 
the  positive  real  axis  in  the  complex  v-plane.  On  displacing  the 
contour  to  the  right  over  a large  number  of  poles  as  shown  in  Figure 
2, we  have  a finite  sum  of  residues  and  a contour  integral  in  which  the 
Legendre  functions  are  replaced  by  their  asymptotic  representations  to 
order  1/v  with  v>>  1.  In  the  finite  sum  of  residues  the  Bessel  function 
is  replaced  by  an  asymptotic  representation,  while  in  the  contour  integral 
its  Sommerfeld  integral  representation  is  used.  The  Bessel  function 
contour  (noted  as  Cy  in  Figure  3 is  displaced  so  as  to  pass  through 
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the  two  saddle  points  at  y = 0 and  y = n in  the  complex  y-plane.  As 
noted  by  Pridmore-Brown,  the  contribution  from  the  poles  that  are  crossed 
in  this  process  and  the  encountered  branch  cuts  yield  the  optical  and 
transition  fields.  The  saddle  point  integrations  give  terms  proportional 
to  (ka)"*l/2.  These  terms  and  the  finite  residue  sum  are  then  transformed 
into  an  equivalent  modified  residue  series.  The  number  of  terms  in 
the  residue  series  is  then  allowed  to  increase  without  limit.  The 
resulting  series  is  summed  in  the  Cesaro  sense  and  represents  the 
field  diffracted  by  the  tip. 

It  should  be  noted  that  in  the  above  method  of  evaluation  of 
these  integrals  (Pridmore-Brown,  1972)  the  condition  ka  >>N>>  1 is 
imposed.  Here  ka  is  the  location  of  the  slot  relative  to  the  tip  of 
the  cone  and  N is  the  number  of  poles  over  which  the  hairpin  contour 
is  shifted  to  the  right.  Subsequently,  N is  allowed  to  increase  without 
limit,  which  requires  that  ka  also  increase  without  limit  if  the 
initial  condition  ka  >>N>>  1 is  to  be  satisfied.  In  fact  it  must  increase 
proportional  to  N^.  Thus,  performing  the  limiting  operation  would 
appear  to  make  all  terms  but  the  optical  term  negligible.  Therefore 
it  is  recognized  that  interpretation  of  the  results  obtained  here  is 
not  clear.  Even  so,  Pridmore-Brown  obtained  numerical  results  for  the 
6-component  of  the  field  that  were  in  very  good  agreement  with  results 
computed  from  the  modal  series  even  for  moderate  values  of  ka.  The 
exact  reason  for  this  good  agreement  is  not  yet  understood,  but  should 
be  investigated  further.  In  addition,  the  possibility  of  using  a 
similar  approach,  but  without  letting  N approach  infinity  should  be 
investigated  in  an  effort  to  extend  the  applicability  of  the  results 
to  even  smaller  values  of  ka. 
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In  this  report,  the  analysis  of  Pridmore-Rrown  is  reviewed,  some 
typographical  errors  are  corrected  and  some  additional  details  are 
supplied.  We  shall  begin  by  first  giving  the  formal  defining  equations 
of  both  components  of  the  electric  field  in  Section  2.  We  shall 
also  point  out  here  the  necessary  corrections  and  omissions  in  the  supporting 
equations  leading  to  the  final  expression  of  the  6-component  of  the 
electric  field. 

In  Section  3,  the  analysis  is  extended  to  apply  to  the  ^-component 
of  the  electric  field,  resulting  in  the  separation  of  this  field 
component  into  transition  and  diffraction  fields.  It  is  found  that, 
unlike  the  6-component,  the  ^-component  is  not  represented  by  an 
optical  term,  whereas  the  diffraction  coefficients  have  the  same 
variation  in  ka  for  both  components.  Throughout  the  above  two  sections 
we  follow  very  closely  the  formulation  and  analysis  presented  by 
Pridmore-Brown  (1972)  and  the  companion  paper  (Pridmore-Brown , 1973). 

The  same  notation  will  be  used  throughout  unless  specifically  stated 
otherwise.  Some  algebraic  details  and  transformations  associated  with 
the  various  integrations  will  be  included  in  the  appendices. 

Also  during  this  program,  the  radiation  patterns  of  circumferential 
and  radial  slot  elements  6.22  wavelengths  from  the  tip  of  a cone  were 
measured  at  a frequency  of  10.38  GHz.  A number  of  the  measured 
patterns  are  shown  together  with  corresponding  computed  patterns 
in  Section  A.  The  computed  patterns  were  obtained  from  the  modal 
series  computer  program.  Except  for  a few  discrepancies  which  are 
explained,  the  agreement  between  measured  and  computed  patterns  is 
very  good  for  both  the  circumferentially  and  the  radially  directed 
element. 
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Figure  1.  Conical  Geometry  and  Slot  Coordinates 


Imv 


Figure  2 


. Displaced  Hairpin  Contour  cv< 
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2.0  GENERAL  FORMULATION  OF  THE  ASYMPTOTIC  APPROACH 

The  geometry  under  consideration  is  a perfectly  conducting  conical 
surface  coinciding  with  the  coordinate  surface  0 = 0Q  of  a spherical  coor- 
dinate system  r ,0,j5  as  shown  in  Fig.  1.  A narrow  azimuthal  slot  located 

at  r = a is  excited  by  an  alternating  voltage 

V0  cosm«5  exp(-io)t)  (1) 

The  external  field  generated  by  this  excitation  is  given  by 

E = curl  curl  (rl^)  + k curl  (rT2)  (2) 

where  and  >2  are  scalar  potential  functions  and  are  given  by  equations 
(3)  and  (4)  of  Pridmore-Brown  (1972).  The  field  components  are  obtained 
from  equation  (2),  giving 

Vo  k 

Eq  - (2iTka)  cosing  exp  [ i (kr  - ir/4)  ] PQ  (3) 

where 

P _ p im2 

0 J0  ~ sin0  sin0n  P20 


JL_ 

2Tri 


vdv 

v2-ii 


A '(0,0  ) i“v 


iv 


J,,(ka) 


A2v(e,eo)  i 


V(ka) 


p 111  (COS0) 

Alv(0,e„)  . -!=U 

PA  (COSeo) 


(3b) 


(3c) 


(4) 


7 


(5) 


P m (cos0) 

a.  (e,e0)  = 3 — 7 

2v  (3/30)  P m (cos0) 

\>-h 


Vx)  = 


(il)*5 


(x) 


and 


E = — / 2nka  m sinmii  exp (i (kr  - n/4) ] P/ 

4>  r p 


where 


P<5  sin0  ?l*  + sin0Q  ^ U 


(6) 


(7) 


(7a) 


vdv  A (0,0Q)  i V Jy(ka)  (7b) 

v2  -l4 

vdv  A ' (0,0O)  iV  ?'(ka)  (7c) 

v2  - 2v  v 

The  integrals  are  taken  over  a hairpin  contour  enclosing  the  positive 
real  axis  of  v > ^ (Primes  denote  derivatives  with  respect  to  the  argu- 
ment or  with  respect  to  the  first  argument  in  the  case  of  two  arguments. 

For  clarification,  we  have  also  added  the  subscript  0,^  to  denote 
quantities  corresponding  to  and  E^,  respectively.) 

The  above  expressions  may  also  be  obtained  by  inspection  from  the 
corresponding  azimuthal  terms  in  the  modal  series  expressions,  noting  the 
exp(jwt)  time  dependence  which  has  been  assumed  in  our  previous  reports. 
This  variation  is  in  contrast  to  the  exp  (-iu>t)  used  by  Pridmore-Brown. 
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This  difference  is  accounted  for  by  taking  the  complex  conjugate  of  the 
modal  series  in  our  previous  reports.  In  addition,  the  slot  dimension  in 
the  previous  reports  introduces  a constant  multiplier  to  the  fields.  Both 
of  these  points  must  be  kept  in  mind  in  any  comparison  between  our  previous 
modal  series  results  and  the  results  from  the  present  analysis. 


2.1  On  the  Transition  Field  of  E 

q 


In  this  section  we  shall  make  the  necessary  corrections  to  the  different 
expressions  given  by  Pridmore-Brown  (1972)  and  in  addition  derive  some 
missing  quantities  necessary  for  a useful  expression  of  the  0-component.  In 
the  following  corrections  the  numbers  of  the  equations  used  by  Pi idmore-Brown 
(1972)  are  retained  for  convenience.  Equation  (13)  should  read 


Bln(e’9o)  = lim(v-vn)  Alv(0,0o) 
v->v 

n 


-m 

p (cos0) 

V~% 


(3/3v)  P (cos0  ) 
v-H  ° 


V = V 


n 


(13) 


Similarly,  the  residue  of  A2n  at  its  nth  pole  is 


B2n(0’V  = llm(^n)  A2M(9’eo> 
P+Pn 


p m (cosG) 

v~h 

(3  2 /a  «oa  u ) p^_^  (cos0o) 

P = V 

n 


(13a) 


The  expression  for  4>n ( y , 0)  below  equation  (23)  should  be 


%(Y,0)  - r"  + 1 *ln  + r“  X2„  + r-"  x3n  + r-"-1  x4n  (23,) 

Equation  (26)  should  be  multiplied  by  (-1) . 
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The  expression  for  the  complementary  error  function  below  equation  (27) 
should  be 


erf c(z) 


exp(~t3)  dt 


with  asymptotic  form 


erfc(z) 


exp(-z^) 


+ 0(z“3) 


(27a) 


(27b) 


Equation  (35)  should  be  multiplied  by  (-1).  In  equation)  (30)  the  arguments 
of  the  f and  g functions  should  be  |u|. 

In  the  evaluation  of  the  first-order  contour  integral,  the  deformed  Cy 
contour  must  also  go  around  the  branch  cuts.  The  branch  cut  contribution 
is  of  comparable  magnitude  to  the  rest  of  the  transition  field  and  must  be 
retained  in  the  final  expression  of  the  field.  A closed  form  expression 
for  the  branch  cut  contribution  is  given  by  the  author  in  the  second  paper 
(Pridmore-Brown,  1973).  The  first-order  transition  field  which  is  evaluated 
for  small  a now  becomes 

T1TF  = /dY  GXp(iRa  C°S  ) ln  ~~  + 3BR}  («>  H («) 


C~fr")1/2  exp  (ka  - J.)  j 


U(U)  + IBR)  (a)H  (ot)  <37a) 


where  H (a)  is  a unit  step  function, 


U(u) 


/T 


/ 


dx  exp(-x^)  ln(w-x)  -lnw 
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and  I (a)  represents  the  contribution  from  the  branch  cut.  The  details 
BR 

of  the  above  evaluation  are  included  in  Appendix  A along  with  further 
description  of  the  U(u)  function  in  terms  of  the  plasma  dispersion  function. 

(1)' 

The  leading  terms  in  the  expansion  of  the  function  A (0,8o)  of 
equation  (20)  in  traveling  waves  are 


(1)’ 
1 


(«’V 


i(*  + *) 

e ° 


(38a) 


where 

n 3 + 4m?  cot  q + A — cot  q 

*8  8 ° 

The  above  result  of  the  expansion  is  used  to  obtain  the  complete 
first-order  transition  field  which  includes  the  contributions  from  the 
branch  cuts  corresponding  to  ot^g  and  a20'  The  first-order  transition  field 
is  then 


2.2  On  the  Diffracted  Field  of  E 

***'  1 1 

In  the  approach  used  by  Pridmore-Brown  the  diffracted  field  is  given 
by  the  Cesaro  sum  (Wittaker  and  Watson,  1944)  of  the  residue  series.  In 
calculating  the  Cesaro  sum,  the  terms  in  the  residue  series  are  modified 
by  adding  and  subtracting  terms  whose  Cesaro  sura  can  be  calculated  in 
closed  form.  This  procedure  results  in  a closed  form  sum  plus  a convergent 
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W 


q 

u 


U 


series  that  can  be  summed  numerically.  The  terms  to  be  added  and  sub- 

» 

tracted  are  determined  by  the  following  method.  The  function  B^n(0,9o) 

is  first  expanded  into  zeroth-order  and  first-order  terms  in  1—  . The 

n 

details  of  this  expansion  are  shown  in  Appendix  B.  The  Cesaro  sum  of 
these  zeroth-order  terms  is  then 


Q(A.0o)  =EQn<e’eo>=£ 


J_  / Sln9o\ 

eo  \sine  / 


(0) 

cos  un  (0Q-0)e 


(C,l) 


-i  /slnQo\^ 
40Q  Uin0  J 


lpa 


ipB 


i 
e 

T 

sina  sing 


(47) 


where 


a = 


B = 


— (tr  -0  +9) 

20o 


20 


(tt+0„  - 0) 


-(m+  H) 


Similarly,  a closed  form  for  the  sum  of  the  first-order  terms  may  be  found. 


R(0,0n)  = L Rn(e.Q0)- 


(47a) 


The  individual  terms,  Rn(0,0o),  in  the  above  series  are  obtained  from  the 

(0) 

original  first-order  expansion  by  replacing  the  magnitude  factor  l/vn 
by  0o/niT,  giving 


Rn(®»®0^  = nTr 


1 /sin9oV2  r 

ITT  ysin0  J nl  s 


sin  v^(0  -0)  -in-,  cos  v^(0  -0) 
no  z n ° 


, , (0)  , 

exp(-i  vn  tt) 


(47b) 


Substituting  (47b)  in  (47a)  and  carrying  out  the  indicated  summation,  we 
find 
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} 
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R(M0) 


where 


{nl(El(’>0'8l  (S-5,)  + n2(gi(a.»4Si(6,oU 


(47c) 


Sl(*,Z)  * -e"  * [in (2  |sina|)  -i(a  . - 3L  )] 

modTr  2 

8l(e,{)  - -e'165  [l„(2|slnB|)  -K^,-  i )] 

5 “ m - \ 

a, 8 defined  below  Equation  (47)  above. 

A complete  expression  of  Ef)  radiation  field  in  tenns  of  the  optical, 
transition,  and  diffraction  fields  may  now  be  written  from  the  above  results 
and  the  results  of  Pr idmore-Brown  (1972,1973).  Thus, 

* I /-COSQ  \ % 

c0  = cos  m4>  j I j [exp(ikR)  + T exp(ikR’)  H(tt  - 0 -0)  +tf] 

( \ sin0  / ° 

exp[ik(r+a)]  , ) 

+ — i — ~L±~  a (0,0  )> 

(ka)^  ° > 


where 


^0  - Efl  -L  [W  COt^  - 0n)l  h 


r-a  cos  (0  - 0 ) 


R'  = r-a  cos (0  + 0 ) 


H(a)  is  a step  function,  (TT(a)  - 1 for  a > 0),^  0. 


1 lkr  ( O)  (1)  / \U 

TF.-e  {irn.iBR(a10).H(a10)-1%iBR  (^i)elkac(Ui0)} 

In  the  above  equation  the  symbols  have  the  following  meanings: 
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r - i (— l ) 

3 + 4m2 


1 “ 4m2 

q cot  0 i-  r—  cot  0 
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(1) 

xBR  M - *= 


ika  cosa[l+H  tan2 


a]  /-  2tt1  \h  (f 

\ ka  cosa  / 


ika  cosa  tan  n 


T (u)  - 

• -u2  0 / 

i erfc  [ 

l-±  r-  | , > 

T~  u 

1 + i 

L \ 

2 / 

TT  U 

sgn(u) 


'10 


l10 


(¥)'  . 

(“)'  • 
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IT  - 0 - 0 

O 


a20 


U(u)  - 


TT  - 0Q  + 0 


CO 

±[ 

/tf  J 


(lx  exp(-x  ) Ln(w-x)  - Lnw 


/lka\  h / ±tt  \ H 

V 2 / “ " ( 2 ) “ 


o_(6,0  ) = (2Tricot0  ) 

m'  o o 


OO  00 

[Q<9,8o)  +E(e'V  Sm+  .ln80  2 s2„l 

L n= 1 ° n=l  J 


’In 


FT  Bln  exp(-ivn7T)  - Qn(0,eo)  - ^(9,0  ) 


v * - % 
n 


Q>  Qn»  R.  Rn  are  given  by  equations  (47),  (47a),  and  (47b), 


Bln'<e,9o> 


-m 


<3%)  Pyn  - 'i  (c°se) 


<5Tr)  p“"  u (=086  ) 

n v - \ u 
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B2n(9’eo> 


-m 

pu  (cosG) 
n - H 


(~r)  p'm 


a jfl'1  P l (cosG  ) 

3p  3 0 Un  ~ ^ ° 


n o 


-m 


-m 


vn » are  the  roots  of  p . (cosfl  ) = 0 and  (3/3fl  ) P , (cosG„), 

v - *5  ° v - h ° 

respectively,  as  in  the  modal  series,  and 


(0) 


= (nn  - tt / 4 + m rr / 2 ) / G 


n = 1,2, 


The  notation  Ln  denotes  the  principal  value  of  the  logarithm  (|arg|  < tt)  , 

It  is  seen  that  the  diffraction  coefficients  am(0,0o)  are  explicit 
fur.ctiqns  of  the  coordinates  0,0q.  It  is  this  property  of  the  diffraction 
coefficients  which  makes  this  analysis  attractive  for  computation  of 
patterns  from  slots  far  away  from  the  tip  of  the  cone. 

In  numerical  calculations  the  seemingly  singular  point  at  a = n/2 
of  the  branch  cut  function,  Ig^  (a),  is  overcome  by  using  the  asymptotic 
form  of  the  defining  expression.  This  gives 


2ne 


ika  cosa 


ka  sina 


which  at  n/2  becomes  2ir/ka  end  can  be  neglected  when  compared  to  other  items. 
Also,  the  function  T(u)  may  be  alternatively  computed  by  the  expression 

T(u)  = [f(|u|)+a(|u|)  + i [f(|u|)-g(|u|)]  - -l±i-  ] sgn(u) 

TT  |u| 

where  f (|u  | ) andg(|u  |)  are  functions  related  to  the  Fresnel  integrals  and  have 
simple  rational  approximations  [Abramowitz  and  Stegun,  1964.] 
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As  it  was  pointed  out  in  the  introduction,  the  above  expression  of 
the  diffracted  field  has  been  obtained  by  allowing  N to  go  to  oo . 

Pdiff  = N - JPlres  + PlsP)' 

This  limiting  process  disregrrds  the  condition  ka  >>  N >>  1 imposed  at 
the  outset,  unless  ka  is  also  allowed  to  go  to  m faster  than  N.  The 
field  expressions  derived  by  using  this  limiting  process  are  then  used 
for  relatively  small  values  of  ka,  ka  > 10,  with  satisfactory  results 
over  most  of  the  range  of  0 (Pridmore-Brown , 1972). 

At  this  point  an  additional  comment  on  the  procedure  of  Pridmore- 
Brown  is  in  order.  In  evaluating  the  zero  order  integral  by  the  method 
of  steepest  descents  he  separates  the  result  into  a so-called  saddle 
point  contribution  and  a transition  field  that  results  when  a pole  exists 
near  the  saddle  point.  The  two  terms  together  remain  finite  when  the 
pole  approaches  the  saddle  point.  However,  separately,  each  has  a singu- 
larity at  the  saddle  point.  In  performing  his  Cesaro  summation,  Pridmore- 
Brown  argues  that  the  series  containing  the  contributions  from  the  saddle 
points  gives  no  contribution  to  the  sum  and  it  is  neglected.  The  remaining 
transition  field  now  has  a singularity  when  the  pole  is  at  the  saddle 
point.  However,  the  Cesaro  sum  of  the  residue  series  provides  the  proper 
singularity  to  cancel  the  singularity  in  the  transition  fields  and  the 
resulting  fields  remain  finite.  A similar  type  of  thing  happens  with 
the  first  order  terms  as  well.  Thus,  though  his  diffraction  fields  and 
transition  fields  each  have  singularities,  the  total  fields  determined 
by  his  method  remain  finite. 
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It  may  be  possible,  however,  to  separate  the  diffraction  field  from 


the  optical  and  the  transition  fields  without  the  need  to  let  N to  to  « 

by  simply  adding  the  non-singular  terms  of  to  the  two  finite  residue 

sums.  The  terms  of  P containing  the  singularities,  (n  = 0 terms),  are 

lsp 

grouped  with  the  transition  field.  In  this  manner  it  may  be  possible  to 
apply  the  technique  to  smaller  values  of  ka  than  would  otherwise  be 
possible.  At  the  same  time  the  singularities  in  the  diffracted  field 
introduced  by  the  limiting  process  will  no  longer  be  present.  This 
alternate  approach  should  oe  examined  further.  Also,  numerical  computations 
should  be  undertaken  to  verify  the  study  results. 


3.0  REPRESENTATIONS  OF  THE  E -COMPONENT 

p 

The  formal  expressions  of  the  ^-component  of  the  electric  field  are 
given  by  equation  (7)  where  the  integrals  are  taken  over  the  hairpin  con- 
tour enclosing  the  positive  real  axis  of  v,  as  was  done  for  the  T -com- 

w 

ponent.  In  comparing  the  form  of  the  integral  expressions  of  this  com- 
ponent with  that  of  E^,  we  see  that  the  only  difference  Is  the  interchange 
of  the  differentiation  with  respect  to  8 of  A,  and  A„  functions.  In 
view  of  this  similarity,  we  may  proceed  witli  an  analysis  completely 
analogous  to  F.  as  was  presented  by  Pr idmore-Rrown  (1972)  and  modified  in 
Section  2.  For  brevity,  however,  we  shall  only  write  down  the  important 
end  results  and  definitions,  keeping  the  same  basic  notation.  The  nota- 
tion and  meaning  of  the  various  quantities  shall  be  the  same  as  for  the 
Eg-component , except  for  the  identifying  subscript  i (for  ^-component)  in 
areas  where  the  distinction  is  necessary  for  clarity. 


3.1  Representation  In  Terms  of  a Shifted  Contour 


As  before,  we  start  by  shifting  the  hairpin  contour  to  the  right  over 

N poles,  where  ka>>N>>l,  resulting  in  a finite  sum  of  N residues  and  an 

integral  over  the  contour  C shown  in  Figure  2.  In  these  contour  integrals 

we  expand  the  angularly  dependent  functions  in  powers  of  retaining  terms 

v 

up  to  0(v  in  the  integrand  as  was  done  for  the  E^-component . 
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Thus 


Pl«*ci 


P 

2<ici 


1 

TiT 


1 

2iri 


11m 

v -»  “> 
o 


11m 

V -►  o» 

o 


J^(ka)  cv  (3.1) 


j'(ka)  dv  (3.2) 

v 


where  the  expansions  of  A (0,9  ) and  A'  (0,0  ) are  the  expressions  in 

■tv  o 2v  o 


the  brackets  and 


J'  (ka) 

v 


kl  + 


J (ka) 

v 

2ka  " 


- J 


v+1 


(ka) 


(3.3) 


From  the  above  it  is  evident  that  the  will  result  in  a first-order 

contour  integral.  P^^,  on  the  other  hand,  will  result  in  a zeroth-order 
contour  integral  arising  from  the  first  term  of  (3.3)  and  a first-order 
contour  integral  from  the  last  two  terms  of  (3.3).  The  first  two  terms 
of  (3.3),  however,  vary  as  ~ as  compared  with  the  other  contour  integrals 
ana  hence  may  be  neglected.  This  is  also  the  reason  why  the  second  con- 
tour integral  in  the  EQ  - component  goes  to  zero  and  does  not  have  any 
significant  contribution  to  the  transition  field.  We  shall  consider  this 
next  order  in  ka  for  the  E^-component  only  so  as  to  gain  some  insight  into  the 

significance  of  each  of  the  terms  in  the  final  E,  expression. 

v 


3 . 2 Zeroth-Order  Contour  Integral 


By  expanding 


into  traveling  waves  valid  for  each  of  the  two  sec- 


cos^0 

tions  of  the  contour  (Imo>  0,  lmv<0)  and  substituting  the  Sommerfeld 
integral  representation  of  the  Bessel  function  we  find  for  the  zero-order 
contour  integral 
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p(0) 

2<*ci 


00 


J (-l)n  <My,6)dv  (3.4) 


where  4 (y,Q)  is  given  by  Equation  (23a)  of  Section  2.1.  For  convenience, 
n 

we  show  in  Figure  3 the  C contour  and  the  deformed  contour  passing 
through  the  saddle  points.  Figure  3 also  shows  the  location  of  the  simple 

po1es  (“iq  » a20  * 01  an  ^ t*ie  zeroth-order  contour  integral  that  must  be  crossed. 

The  indicated  ’'ranch  cuts  are  associated  with  the  singularities  (a  . a ) 

10  ?0 

of  the  first-order  contour  integrals  to  be  considered  in  the  next  section 
and  are  taken  parallel  to  the  imaginary  axis. 


We  may  now  write  the  field  arising  from  crossing  the  poles  and  the  zeroth- 
order  transition  field  from  the  integration  along  the  deformed  contour 
passing  through  the  saddle  points  as 


, (0) 
2$c 


i ' Ski  (Si2)  i «-lka"S(e°-9>+  r e-iKacos(V9)„(,.ve) 


elkaI'T(u10) 


} 


(3.5) 


where  T(Uj-)  is  defined  as  before.  As  in  the  E case,  the  contributions 

from  the  saddle  point  at  y = it  for  and  u^^  cancel  one  another  since 

T (u)  - -T  (-u).  It  is  interesting  to  note  that  these  fields  differ  by 

the  factor  — from  the  corresponding  fields  of  the  E -component, 
ka  9 


20 


w 


fi 


I 1 

m 

* 

I I 

I , I 


We  have  purposely  considered  the  terms  proportional  to  (ka)  so  that 
we  may  compare  the  different  types  of  fields  with  the  corresponding 
fields  of  the  G-component.  It  is  immediately  evident  that  the  optical 
field  in  this  case  is  not  at  all  significant.  Therefore,  these  terms 
must  be  neglected  if  the  9-component  and  the  ^-component  are  to  be 

-h 

computed  to  the  same  order  in  ka , i.e.  (ka) 


3.3  First-Order  Contour  Integral 


In  the  evaluation  of  the  first-order  contour  integrals  again  we 


make  use  of  the  traveling  wave  expansion  of  the  function  — 1- — and 

sinV0 

^ for  the  shifted  hairpin  contour  Cy 


In  taking  the  deformed  C 


cos^0  ' v Y 

contour  through  the  saddle  points  we  encounter  branch  points  associated 

with  the  logarithmic  singularities  of  the  exponential  integrals.  The 

contour  integrals  associated  with  (ka)  in  and  give 

contributions  identical  to  those  of  the  E0-component  and  can  be  immed- 

. (1) 

lately  written  by  analogy.  However,  the  contour  integral  of  P2<j>ci 
arising  from  the  second  factor  of  (3.3)  varies  as  1/ka  and  need  not  be 
considered  further.  Thus, 


(1 

l$c 


> i /sln"o  \ 
i ' 4,2  (Sl"«  / 


(1)  (O 

-r  x-  <aio>  + IBR  (a20) 


BF. 
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i(ka  - tt/4) 


ka 


U(U10> 


(3.6) 


(D 


2(Jici 


a 0 


(3.7) 
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which  is  significant 


The  only  first-order  contour  integral  of  P2<j>ci 
(magnitude  proportional  to  (ka)_>s)  is  the  one  arising  from  the  last 
term  of  (3.3).  This  contour  integral  has  the  form 

ikacosy  iy 


( 1 ) B 
[TF 


/ 

Cy 


dy  e 


^ (ivQ  (a  - y))-  E^IVqO) 


(3.8) 


differing  from  the  previous  first-order  contour  integrals  by  the  factor 
eiY  in  the  integrand.  The  same  process  of  integration  through  the  saddle 
points  at  y = 0,  y « tt,  and  around  the  branch  points  is  carried  out  as 
before,  giving 


(1) 

) 

2<}>c 


B = -1  / slneoVi  [j 

1 mr2  \sine  J L 


(1) 


riBR^(a10>H(a10>  + 


W<4«>  -rJE 'Uka-VMW 


rr—i  -i (ka-ir/i*)  r * * 1 

■J  — e U.  (0  ) - U (u  ) 

" ka  L * 20  ^ 30  J 


(3.9) 


It  is  pointed  out  that  in  the  above  expression  we  have  deleted  the  phase 
i 


factor  e 
.0) 


2ka  . The  functions  U,  (u)  and  the  branch  cut  contributions 


I (a)  differ  from  their  counterparts  in  the  E0-component  but  have 
essentially  the  same  properties.  Unlike  the  case  of  the  Eg-coinponent , 
the  above  contour  integrals  through  the  saddle  point  at  y = n for  u 

and  u result  in  a net  contribution  to  the  field  since  U.  (u)  is  not 

3°  * 

an  even  function  of  u and  hence  U,  (n  ) ^ Ilf  (C  ) . If  upon  further 

v 20  <P  30 

examination  this  contribution  is  found  to  be  insignificant,  then  these 
factors  will  not  be  included  in  the  computation  of  the  field. 
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3. A Evaluation  of  the  E .-Diffracted  Field 

The  diffraction  field  of  the  E -component  is  given  by  the  sum  of 
the  residue  series, 


-IP,  + i P 
s inG  * s in© 


2i|>res 


4>diff 


(3.10) 


o 

The  residue  series  P . and  P0,  are  found  from  the  expansion  of 

lvres  2<f>res  1 

the  functions  B^(G,  0Q)  and  ®o^ * resPectively , and  the  asymptotic 

representation  of  the  Bessel  functions  J (ka)  , Jv+1(ka).  Thus, 
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(3.11) 


(3.12) 


where  Bln  (Q,Go)  and  b2n(e,eo^  arG  8iven  by  (13)  and  (13a)  in  Section  2.1. 

Both  residue  series  converge  in  the  ordinary  sense.  This  summation  is 

accomplished  by  expanding  the  coefficients  B.  (6,  G ) in  powers  of  l/v 

in  ° n 

f ! _ 

and  ®2n  (®’Go>  bn  Powers  °f  (0)  and  retaining  the  leading  term  in  each 

Wn 

expansion.  However,  on  retaining  only  terms  which  vary  as  (ka)  *,  as 
in  the  Ey-component , we  are  left  only  with  two  convergent  first-order 
series  denoted  by  11^(6, Gq)  and  R2^  (G  ,©0)  and  defined  below: 
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where 


and 
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(3.15) 


(3.16) 


are  the  individual  terms  in  the  two  series.  The  details  of  the  expansion 
cf  Bln(0,Oo)  and  B2n(0,0o)  leading  to  Rj^ , R1(J) , R2n^ » R2$  are  included 
in  Appendix  B.  The  diffraction  field  may  be  now  found  by  the  sum  of  the 
two  series  as 
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and 


Rln^  (0*eo)  * R2n<}i^e,eo^  are  defined  by  (3.15)  and  (3.16). 
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We  also  have: 


= (n 


+ -)  - 

2 «o 
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(3.21) 


The  functions  g^oi.C)  , g1(B,I)  are  obtained  from  g^a.O, 
g^ (6*0  by  replacing  r with  (C+l) . All  other  quantities  have  the  same 
meaning  as  for  the  Eg-component  and  have  been  defined  in  Section  2.3. 

We  are  now  in  a position  to  express  the  radiation  field  of  the  d> -component 
as  a sum  of  transition  and  diffraction  fields  as  was  done  for  the 
Eg-component . In  summing  the  various  fields  and  normalizing  to  the  same 
factors  as  for  the  Eg-component,  we  find 
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Also,  R , R.  , R„  , R_  are  defined  by  equations  (3 . 13)  ^through  (3.16). 
10  in0  /0 

U (u)  and  I„D  (a)  are  given  in  Appendix  A.  All  other  quantities  have  the 

$ dK(Ji 

same  meaning  as  for  the  Eg  and  may  be  found  in  Section  2.  As  for  the 
Eg-component , for  a near  ir/2  the  asymptotic  form  of  (A-12)  must  be  used  for 


(1) 

BR$  1 

i (ka  cos  a + a ) 


the  numerical  computation  of  Tgj^ot).  This  is  easily  shown  to  be 


I(1)(a)  ^LS 

BR$  ka  sin  a -1 
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Since  this  is  of  order  0(l/ka),  the  contribution  of  the  branch  cut  in  the 
vicinity  of  a « tt/2  can  be  neglected  entirely  without  affecting  the 
results  significantly.  It  should  be  noted  that  the  two  diffraction 
coefficients  of  the  0 and  (5— components  have  the  same  variations  in  ka 
but  in  general  will  have  different  values  at  a given  point  in  space. 
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4.0  MEASURED  AND  COMPUTED  ELEMENT  PATTERNS  ON  A CONE 


The  computer  program  which  has  been  developed  for  conical  array 

pattern  analysis  has  been  modified  and  executed  on  the  Control  Data 

Corporation  computer  (CDC  6600).  The  CDC  6600  computer  is  preferred 

because  of  its  greater  accuracy  and  higher  speed  over  the  IBM  370/158 

computer  previously  used  in  the  pattern  analysis  problem. 

During  the  course  of  the  development  of  the  modal  series  program 

the  accurate  computation  of  the  Associated  Legendre  function  P^  (cosG) 

was  questionable  over  certain  ranges  of  0.  Specifically,  the  chosen 

TT  IT 

expansion  was  valid  for  — < 0 < 5 — but  it  was  used  for  computation  of 

6 6 

9 

the  function  for  < 0 <tt  by  first  computing  Pm  (cos0)  with 

7 2+2m  v 

v ~ 72+2m.  Backward  recursion  was  used  to  obtain  the  value  of  the  function 
of  the  desired  degree.  Some  details  in  checking  the  accuracy  of  the 
results  in  both  the  upper  and  lower  extension  limits  of  applicability 
of  the  expansion  are  included  in  the  first  Quarterly  Report 
(Bargeliotea , 1974). 

The  computed  patterns,  using  the  first  fourteen  modes  of  the  modal 

series,  have  been  normalized  to  the  largest  value  computed  for  both 

polarizations  for  either  the  circumferential  or  radial  case.  Computed 

and  measured  elevation  cuts  are  shown  on  an  extended  0 -scale  to  give 

a single  view  of  the  elevation  cuts.  Since  the  cone  under  study  has  a 

o 

10°  half-angle,  the  position  of  0*80  was  selected  for  the  computed 
azimuthal  cuts.  This  corresponds  to  the  broadside  of  the  conical  surface. 
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Measurements  were  made  on  an  experimental  10°  half-angle  cone  at 
a frequency  of  10.38  GHz  with  the  slots  6.22  wavelengths  from  the  cone 
tip.  This  slot  location  corresponds  to  ka  = 39  radians  from  the  cone 
tip  for  which  a complete  set  of  patterns  have  been  computed  with  the 
modal  series  program.  In  both  the  circumferential  and  the  radial  case, 
the  slots  were  fed  by  a half  height  X-band  waveguide.  Where  possible, 
computed  pattern  values  are  shown  on  the  same  scale  for  comparison. 

A . 1 Radiation  patterns  of  circumferential  slot. 

Figure  A shows  the  Q-polarized  patterns  for  *5  = 0°  and  <5  = 180° 
together  with  corresponding  computed  patterns.  There  is  excellent 
agreement  between  the  measured  and  computed  patterns  throughout  the 
cone  tip  region  and  the  broadside  region  where  both  patterns  are  quite 
uniform.  The  measured  pattern  drop-off  for  <5  = 0°  and  9 greater  than 
1A0°  is  due  to  shadowing  at  the  base  of  the  cone  by  the  absorbent  material 
in  which  the  cone  was  set.  Similar  drop-off  is  observed  in  the  measured 
pattern  for  - 180°  for  the  same  0 values  except  that  this  is  not  as 
noticeable  since  the  computed  pattern  in  this  region  is  at  a much  lower 
level  relative  to  the  <5  = 0°  pattern.  Measured  Q-polarized  patterns  are 
also  shown  for  = A5°,  <5  = 225°  and  i = 90°,  = 270°  in  Figures  5 and 

6,  respectively.  In  Figure  6 the  corresponding  computed  pattern  values 
are  also  shown.  There  appears  to  be  a 2°  difference  in  the  position  of 
the  null  between  the  computed  and  measured  patterns  but  this  is  probably 
caused  by  inaccurate  alignment  between  the  transmitter  and  receiving 
antennas.  This  also  explains  the  slight  asymmetry  in  the  measured  pattern. 
The  shadowing  by  the  absorbent  material  for  9 greater  than  1A0°  is 
evident  in  all  measured  patterns.  The  above  two  patterns  also  show  that 
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the  pattern  level  at  4 *»  135  is  approximately  10  dB  below  that  of  = 90  . 

The  {^-polarized  patterns  in  elevation  are  shown  in  Figures  7 through 
9 for  the  same  azimuth  positions  as  for  the  9 polarized  patterns.  Azimuthal 
patterns  for  both  polarizations  at  9 = 80°  are  shown  in  Figure  10. 

Excellent  agreement  is  seen  between  computed  and  measured  patterns  of 
Figures  9 and  10  considering  the  fact  that  computed  pattern  values  were 
made  every  5 degrees  in  9 and  10  degrees  in  i. 

4.2  Radiation  patterns  of  radial  slot. 

The  tip  section  of  the  1C°  half-angle  cone  was  modified  to  accommodate 
a radial  slot  6.22  wavelengths  from  the  cone  tip.  The  edges  of  the  wave- 
guide were  carefully  taped  to  obtain  one-half  wavelength  at  the  operating 
frequency.  A photograph  of  the  cone  mode  with  a radial  slot  and  absorbent 
material  surrounding  the  base  of  the  cone  is  shown  in  Fgiure  11. 

The  computed  patterns  shown  together  with  the  measured  patterns 
have  been  normalized  to  the  largest  value  computed  for  both  polarizations 
using  the  first  fourteen  modes  of  the  modal  series.  This  largest  value 
was  computed  at  the  field  point  P(«5,9)  * P(0°,80°).  As  expected,  this 
point  is  broadside  to  the  conical  surface  and  lies  on  the  plane  normal 
to  tht  conical  surface  and  containing  the  radial  axis  of  the  slot  element. 
For  meaningful  comparison  of  measured  and  computed  patterns,  the 
reference  level  cf  the  measured  patterns  was  also  taken  at  broadside  at 

<5  = 0°. 

Figure  i2  shows  the  »5-polarized  patterns  for  6 = 0°  and  $ = 180 
together  with  corresponding  computed  patterns.  There  is  excellent  agree- 
ment between  the  measured  and  computed  patterns  throughout  the  cone  tip 
region  and  the  broadside  region  where  both  patterns  are  quite  broad. 

3-* 
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igure  10.  Measured  and  Computed  (M  = 13)  Patterns  of  A/2 
Circumferential  Slot  for  0 and  d Polarizations  at 


Fig.  11  Test  Model  of  10°  Half -Angle  C.one  with 
Slot  6.  22A  from  Tip 
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Figure  12.  Measured  and  Compu 
of  A/2  Radial  Slot  for 


Measured  ei-polarized  patterns  are  also  shown  for  d 10°,  d 190°, 

® - 40  , d 220  , © 90  and  d 270  in  Figures  13  through  19.  Computed 

pattern  values  are  also  shown  in  the  same  figures  for  comparison.  Figure 
l 5 shows  a 3°  difference  in  the  position  of  the  null  between  the  computed 
and  measured  patterns  but  this  is  probably  caused  by  misalignment  between 
the  transmitter  and  receiving  antennas.  It  would  also  explain  the  slight 
asymmetry  in  the  measured  pattern.  The  shadowing  by  the  absorbent 
material  for  0 greater  than  140  is  again  evident  here.  The  above  figures 
also  show  that  there  is  a difference  of  4 dB  in  the  maxima  of  d 0°  and 
d *0  patterns,  whereas  the  pattern  maxima  at  d 180°are  approximately 
14  dB  below  that  at  d 0 . As  mentioned  earlier,  the  maxima  occur  at 
the  broadside  to  the  conical  surface. 

I he  Q-polarized  patterns  in  elevation  are  shown  in  Figures  16  through 
19  for  the  same  azimuth  positions  as  for  the  ^-polarized  patterns.  These 
figures  clearly  show  that  unlike  the  patterns  of  a circumferential  slot  on 
a similar  cone,  the  cross-polarization  is  well  below  the  dominant  pola  r- 
ization  at  the  same  azimuthal  angles.  The  shadowing  by  the  absorbent 
material  at  the  base  of  the  cone  is  clearly  identified  in  these  patterns. 

Since  the  pattern  levels  of  this  polarization  are  approximately  18  cl  13 
below  the  reference  level,  the  misalignment  of  the  transmitter  and  receiver 
antennas  is  also  more  apparent. 

Azimuthal  cuts  at  various  angles,  including  the  broadside  of  the  antenna 
element,  are  shown  for  both  polarizations  in  Figures  20  through  24.  From 
Figure  21  it  is  seen  that  the  measured  broadside  pattern  of  the  d-polariza- 
tion  exhibits  a broad  maximum  at  d — 0 and  is  in  very  good  agreement  with 
the  computed  pattern.  As  stated  earlier,  this  has  been  selected  as  the 
reference  point  for  all  patterns.  In  the.  0-polarized  azimuthal  cuts,  where 
the  levels  are  considerably  lower  than  the  reference  level,  the  shift  of 
the  null  between  measured  and  computed  patterns  is  more  apparent.  At 
other  points,  however,  the  agreement  between  measured  and  computed 
values  is  very  good.  The  O-polarized  pattern  levels  are  in  all  cases  at 
least  18  dB  below  the  reference  level. 
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ure  13.  Measured  and  Computed  (M  = 13)  d-Polarized  Patterns 
of  A/2  Radial  Slot  for  6 10°,  6 - 190° 


juauodujo;}  ‘gp  ui  jaMOj 


igure  17.  Measured  9-Polarized  Patterns 

of  A/2  Radial  Slot  for  6 - 10°,  d = 190 


figure  18.  Measured  and  Computed  (M  = 13)  9-Polarized  Patterns 
of  A/2  Radial  Slot  for  d = 40°,  d = 220° 
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igure  19.  Measured  and  Computed  (M  = 13)  ©-Polarized  Patterns 
of  X/2  Radial  Slot  for  d>  = 90°,  d = 270° 


Figure  20.  Measured  and  Computed  (M  = 13)  d-Polarized  Patterns 
of  A/2  Radial  Slot  for  0 = 0°  9=1° 


5.  RECOMMENDATIONS  FOR  FURTHER  INVESTIGATIONS 

Computed  element  patterns  using  the  normal  mode  series  computer 
program  have  shown  excellent  agreement  with  measurements  on  an  ex- 
perimental model  with  both  types  of  slots.  Using  the  model  series 
for  the  computation  of  exact  patterns  the  number  of  modes  required 
for  convergence  of  the  series  increases  as  the  distance  of  the  slot 
location  from  the  tip  of  the  cone  increases.  Computational  difficulties 
from  accumulated  round-off  errors  may  also  affect  the  accuracy  and 
limit  the  applicability  of  the  modal  series  program  to  element  positions 
not  too  far  from  the  cone  tip. 


We  also  have  now  complete  expressions  of  the  Eq  and  E^  radiation 
fields  in  terms  of  optical,  transition  and  diffraction  fields.  Furthermore, 
it  is  seen  that  the  diffraction  coefficients  of  both  components  are  only  a 
function  of  the  angular  coordinates  and  not  the  location  of  the  radiating 
element.  Ihe  diffraction  fields  for  any  element  location  may  then  be 
found  by  the  product  of  these  diffraction  coefficients  and  the  factor  con- 
taining the  information  of  the  location  of  the  radiating  element.  The 
diffract' on  coefficients  need  be  computed  only  once,  thus  eliminating  the 
repetition  of  lengthy  computations.  The  optical  and  transition  portions  of 
the  fields,  on  the  other  hand,  are  given  in  terms  of  functions  that  are  more 
suitable  for  numerical  computation  than  the  modal  series  expressions. 

In  order  to  utilize  effectively  the  results  of  the  above  approximate 
approach  for  computation  of  radiation  patterns,  the  expressions  of  both 
polarizations  must  first  be  programmed  for  use  by  the  digital  computer. 
The  new  program  should  be  built  in  "blocks"  as  much  as  possible,  so  as 
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to  utilize  certain  existing  computing  routines  but  more  importantly  to 
facilitate  the  integration  with  the  modal  series  computer  program.  Such 
an  integration  of  the  two  computer  programs  will  be  necessary  since  the 
approximate  expressions,  by  their  nature,  fail  to  yield  satisfactory  re- 
sults in  the  regions  near  the  cone  axis.  Also,  unlike  the  modal  series, 
the  asymptotic  approach  cannot  be  properly  applied  to  radiating  elements 
near  the  tip  of  the  cone.  Thus,  the  combination  of  the  two  programs  will 
tend  to  complement  one  another. 

Once  the  newly  developed  computer  programer  is  tested  for  general 
accuracy  and  proper  behavior  in  critical  regions,  it  will  provide  the 
means  for  computing  radiation  patterns  of  both  polarizations  from  slots 
on  a cone.  The  testing  should  involve  computation  of  patterns  for  identi- 
cal slot  locations  that  were  computed  by  the  modal  series  and  verified  by 
measurements.  In  addition  to  the  above  numerical  analysis,  the  impli- 
cations of  the  limiting  process  discussed  in  Section  2 should  be  investi- 
gated further. 

The  analysis  technique  outlined  above  would  allow  the  computation 
of  patterns  to  verify  the  equivalence  principle  pattern  synthesis  technique. 
Additionally,  the  technique  could  be  easily  extended  to  analyze  mutual 
coupling  and  impedance  of  radiators  on  a cone  from  a rigorous  approach 
and  provide  a check  on  the  results  obtained  by  the  geometrical  theory  of 
diffraction. 
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APPENDIX  A 


EVALUATION  OF  CONTOUR  INTEGRALS 


In  the  evaluation  of  the  first-order  integrals  we  are  confronted 
with  integrals  of  the  following  form. 


I * f elka  cos  Y In  -3-  dY 
J a-y 


(A-l) 


In  this  expression  is  the  deformed  contour  (Figure  A-l)  passing 
through  the  saddle  point  Y = 0.  In  addition,  when  evaluating  the 
contour  integral  of  Eq.  3.8  for  small  a,  the  following  integral  occurs. 


i k a cos  Y iY  , v 

f e e In  -£L_  dl 

J , a-Y 

C 

Y 


(A-2) 


The  first  integral  may  be  obtained  as  a special  case  of  the  second  so 
that  only  the  second  integral  will  be  discussed  here.  When  a>0  the 
deformed  contour  must  include  a section  along  the  branch  cut.  When 
<*<0  the  branch  cut  lies  entirely  below  the  deformed  contour  and  there  is 
no  integration  along  the  branch  cut.  When  ct>0  the  path  may  be  subdivided 
into  three  parts  as  illustrated  in  Figure  A-2.  The  integral  may  then  be 
written  as  follows: 


i(k  a cos  y+y) 

In  (-2-)  dy 
a-Y' 


(A-3) 


i 


Along  C ? the  logarithm  of  a-y  differs  by  -12*  from  the  value  It  would  have 
have  if  there  were  no  branch  cut.  Consequently,  the  integrals  along  C j 
and  C^2  can  be  treated  as  an  integral  taken  along  both  curves  as  though 
there  were  no  branch  cut  plus  an  integral  along  C^2 . The  deformed  contour 
integral  may  now  be  written  in  the  following  form. 


i (k  a cos  y + y) 


In  (— )dY 
a-y 


f i (k  a cos  y + y) 

/ e Ln  (—2—)  dy 

t a-y 

C C 
Y1  + Y2 


(A-4) 
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i(k  a cos  y+  y) 

? dy  + 


i(k  a cos  yd  y) 
i ln(_2_)dy 

a-y 


where  arg  Ln(a-y)  is  restricted  between  -if  and  it  . 

The  integral  over  the  branch  cut  is  obtained  as  follows: 

a y 

/i  (k  a cos  y + y)  a | i(ka  cos  y + y ) 

e Ln  (^)dV  + I e Ln(^ 

Yl  ,Y,  A 

fi 

I i(k  a cos  y + y ) 

+ 12tt  e dY 

J a 


(A-5a) 


Y 

/'  i(ka  cos  y + Y ) 
e d' 


(A- 5b) 


where  Re  y^*  a 
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The  deformed  contour  may  now  be  written  in  the  following  form: 


t 


i(k  a cos  y + y ) 
e ln(— 2— )dY 

oTy 


/ 


i(k  a cos  y+  Y ) 
e Ln  (_2_J  dY 

“-Y 


C ,+C 
Y1  Y2 


(A-6) 


i(k  a cos  Y+  Y) 
+i2ir  / e dY 


/* 


where  C3  is  a contour  from  a vertically  to  C and  then  along  C as 

Y2  y2 

shown  in  Figure  A-2  . Since,  as  Iiy  ei  ka  cos  y approaches 

zero,  for  Re  y between  0 and  JL  , the  contour  C,  may  be  deformed  to 

2 3 

coincide  with  the  branch  cut  down  to  Im  y->  - Then  the  total  deformed 
contour  integral  has  the  following  form 


l 


l 

1 
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I: 


i (k  a cos  Y + Y) 
e ]n(— )dY 

ot-Y 


I: 


i[k  a cos  y+y] 


Ln(_)dY 

a-Y 


C . + C , 
Yl  Y2 


(A- 7) 
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i(k  a cos  y + y) 
e dY 


Br0 


where  the  subscript  Br^  indicates  a vertical  path  along  the  branch  cut 
from  a to  a - i <*>  . 

rru  u , (D 

The  branch  cut  integral,  denoted  by  1 , may  be  obtained  as 

Bry 

follows. 

Let  a new  variable  y be  defined  by  the  following  relation: 


Y " a -i  y 


(A-8) 
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Then  IBr  takes  on  the  following  form: 


7’  (a)  ‘ 2 


00 
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i k & cosy  iY 


e dy , 


(A-9) 


where  y = (a  - iy) . 

On  using  the  relationships  cos  ( a - iy)  * cos  acosh  y + i sinasinhy 
sinh  y a y 

cosh  y = 1 + y /2  f the  following  expression  results: 


0) 

I (a)  = 2tt 
BrA 


i(k  a cos 


a + a ) f i~-  cos  at  [y2  + i 2 tan  a'  y]  (A-10) 
*T)  dy 


where 


tan  a'  = tan  a - 


k a cos  a 


which  on  completing  the  square  in  the  exponential  becomes 


(1)  i(ka  cos  a + a +~  cos  a tan2  a* ) 

I (u)  - 2ir  e 2 
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The  integral  in  (A-ll)  can  be  put  into  the  form  of  r complementary 
error  function  integral  by  the  following  transformation: 


-t‘ 


i-  cos  ot[  y + i tan  a'  ]1 2 
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Only  the  Integral  over  C + C remains  to  be  evaluated.  Let  y 

Y1  Y 2 

the  real  part  of  y and  let  y"  be  the  imaginary  part  of  y.  Then 
C + C 

Y*  Y2  sinh  yM  = - tan  y' 

To  perform  the  evaluation  the  following  substitution  is  made 

v = sin  y'  / (cos  y')1  2 

Then  on  C + C the  following  relationships  hold 
Y1  Y 2 


ika  cosY  ika  -k  a v2 

e * e e 


be 

on 


(A-l3a) 


(A-13b) 


(A-13c) 


is  very  small  except  when  v is  near  zero. 


If  ka  is  large  then  e 
Then  Eq.  (A-13a  and  b)  reduce  to 


/T  v 


-i  f 

dY  ~ e A /2~dv 


If  the  additional  substitutions 

x 


(A-1V) 


are  used,  the  following  expression  results 


/ 


i (k  a cos  y + Y ) 
e Ln(— — ) dY 


C ,+C 
Yl  Y2 


a-Y 

— i[ka-ib2-  -2-] 


v^F 


/ 


-(x-ib)2 

e Ln  ( w ) dx 

vw-x 


(A-15) 


y^2ka 


where 


Let  U, (u)  be  defined  as  follows: 
<P 


U (u) 
<P 


1 

jr 


(x-ib) 


dx 


(u) 


- (x-ib) 2 

e Ln(w-x)  - Ln  (w) 


(A-16a) 


(A-16b) 


When  a<0  the  deformed  contour  does  not  deviate  about  a branch  cut. 

In  that  case,  i.e,  when  Im  w < 0 only  U^(u)  contributes.  It  is  necessary, 
therefore,  to  determine  methods  for  evaluating  (u) . Because  of 
characteristics  of  the  logarithm,  different  series  expansions  must  be 
used  for  w in  the  upper  half  of  the  complex  plane  and  for  w in  the  lower 
half  of  the  complex  plane.  There  will  be  a discontinuity  as  w crosses 
the  real  axis.  To  determine  the  various  expansions  consider  the  integral 
portion  of  Eq.  (A-16b). 


Fb(w) 


1 


-(x-ib)2 
e Ln 


(w-x)dx 


(A-17) 
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! 


II 

This  function  is  analytic  in  the  upper  half-plane  and  in  the  lower  half- 
plane. The  upper  half-plane  will  be  considered  first.  If  Itn(w)  is 
greater  than  zero  Ffe(w)  is  differentiable  and  it  is  possible  to  write 
the  following  expression 
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This  function  is  analytic  in  the  entire  complex  plane  including  w = 0. 
In  Eq.  (A-18)  the  lower  limit  of  the  integral  may  be  set  equal  to  zero 
and  F^(ie)  may  be  replaced  by  Fb(i0+)  where 

Fb  (i0+)  = lim  Fb  (it)  (A-23) 
e-K) 


This  limit  is  given  by  the  following  expression 


Fb  (i0+) 


-(x+ib)2  -(x-ib)2 

(e  +e  ) In  x dx  + i JL  [l+erf(ib)] 

2 


(A-24) 


and  is  obtained  using  the  condition  -*<arg  lnz  < tt.  The  integral  can 
be  evaluated  numerically.  For  such  an  evaluation  an  integration  by 
parts  is  first  performed  to  remove  the  logarithmic  singularity.  The 
resulting  expression  is  as  follows: 

F(i0+)  = 1(b)  + i j [l  + erf  (ib)]  (A-25) 

where 

l r -(x+ib)2  -(x-ib)2 

1(b)  ; JZ.  J [e  +e  ]lnxdx 

0 


_2_ 

/r 


r°  -(x+ib)2  -(x-ib)2 

J [(x+ib)e  +(x-ib)e  ] x(lnx-l)dx 
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On  using  Eq.  (A-25)  and  integrating  Eq.  (A-22)  the  following  expression 
for  Fb(w)  is  obtained  when  w is  in  the  upper  half-plane. 


i 

t 


:» 


I 


1 

II 


\ 


t S 


[] 


E 


Fb(w) 


Kb) 


+ i i [1 
2 


°°  nr  2n+2  2n+2 

(-l)n[(w-ib)  -(ib) 

erf (w-ib) ]+  ^ 1 • 3* 5* • • (2n+l)  (n+1) 

n*0 


]2n 


Im  w>  0 
(A-26) 


If  w is  in  the  lower  half-plane  then  instead  of  Eq.  (A-20)  the  following 
relationship  is  used. 


1 

x-w 


i(x-w)y 
e dy 


Im  w<  0 (A-27) 


If  an  analysis  is  carried  through  in  a manner  similar  to  the  preceding 
one  the  following  expression  results  for  (w)  in  the  lower  half-plane. 


, r . £ <-l)n[v-lb)  2n+1-(lb)2nl'2]  2n 

F (w)  - Kh)  -i  - [l-erf  (w-ib)]  + £ 1 • 3 • 5 • • • (2n+l)  (n+1) 

1 n-0 


Im  w<  0 
(A-28) 


The  function  U (u)  and  U(u)  are  then  given  by  the  following  expressions 


U^u)  ■ Ffe(w)  - lnw 


(A-29) 


U (u)  - lim  UKu) 
b+0 


(A- 30) 
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A 


It  may  be  shown  that  U(u)  is  an  even  function  of  u.  In  addition 


lim  1(b) 


(A- 31) 


where  C is  Euler's  constant. 


C - 0.577215665. 


(A-32) 


When  Im  w >0  the  branch  cut  integral  must  be  added.  When  Im  w < 0 
there  is  no  branch  cut  integral.  Eqs.  (A-26)  and  (A-28)  permit  computation 
of  U^u)  for  reasonably  small  values  of  |w  - ib  | . For  large  values  of 
|w  - ib j an  asymptotic  expansion  is  desired.  Such  an  expansion  can  be 
obtained  directly  from  the  integral  expression.  To  accomplish  this 
expansion  rewrite  Eq.  (A-16b)  in  the  following  form: 

2 

U.  (u)  = — f e-^x_ib^Ln((w-ib)  - (x-ib))  dx  -Ln  w (A-33) 

A" 


As  long  as  w - ib  is  not  pure  imaginary,  then  as  |w-ib|  increases  without 
limit,  the  main  contribution  to  the  contour  integral  is  from  U^u).  In 

turn,  the  main  contribution  of  U, (u)  will  come  from  the  region  where 

’ <P 

Re(x-ib)  is  a minimum,  i.e.,  where  x=Re  b«  |w-ib|.  Then  the 

logarithm  may  be  expanded  in  powers  of  and  the  resulting  expression 

w-ib 

may  be  integrated  term-by-term.  This  process  leads  to  the  following 


expression: 

UKu)  ' Ln  (w-ib)  - Ln  (w)  - X-  ~ — hr 
v n*l  n(w-ib) 


l r -(x-ib)2 

— r=-  I (x-ib)ne  dx 

n /it 


(A-34) 


- Ln  (w-ib)  - Ln(w)  - 

n=l  2n(w-ib)^n2n 


w-ib  -*•  00 


4 


When  b ■ 0 the  following  expansion  for  U(u)  results: 


u(u>  ~ - E 

n“1  2n  w2n2n 


w -*■  °° 


(A-35) 


The  total  integral  may  then  be  written  as  follows: 


ta  cos  y 


+ iy] 


In  dy 

a-y 


i(ka-ib2-  y) 

I,_  («)  H (a)  - JSL  e 4 lU(u) 

BR  1 ka 


For  large  values  of  ka  the  b2  in  the  exponent  may  be  omitted. 


(A-36) 


il 


An  alternative  derivation  of  the  contour  integral  in  which  the 
branch  cut  is  taken  from  a along  the  negative  real  axis  results  in 
a slightly  different  expansion  for  the  integral  valid  for  ct  sufficiently 
small  that  cos  a ~1  - a2/2.  The  expansion  is  as  follows: 


i(ka  cos  y+y) 

; In  (-^—)  dy 

a y 


rr-  i (ka  — ibz-  — ) p , . 

- e 4 jj(b)  - i J [1  - erf  (w-ib)]  (A-37) 

w t 0 

+ Y"  (-l)n  [(w-ib)2n+2-(ib)2n+2]  2n  - In  wj 
" 1*  3 • 5 * • • 


(2n+l)  (n+1) 


This  expression  is  valid  for  Im  (w)  both  positive  and  negative  and  agrees 
with  the  previous  results  as  w-*-0. 
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When  a is  near  the  saddle  point  at  y“tt  then  the  integral  of  interest 


is  as  follows. 


/*i[  ka  cos  Y + y] 

> ln  S dY 


where  C_^  passes  upward  through  the  saddle  point  at  tt.  This  integral 
leads  to  the  following  integrals 


f i[ka  cos  a cosh  y + i sin  a sinh  y]  i(a-  iy) 

1 * 2tt  / e e dy  H(tt-cx) 


-i(ka  - ib 


Me-  e 


>2  r i f 
' * 1 r 


■°°-(x+ib*) 


Ln  (-- — ) dx 
w -x 


(1)*  [T  -i(ka  - ib  --) 

e 0*(.1 


u . ( y > 
4* 


where 


tt  - a 


w = \T 


e (tt  - a) 


APPENDIX  B 

EXPANSION  OF  B'  (0,  0 ),  B.  (0,0  ), 

In  o in  o 


B 


2n 


(0.  eo) 


Expansion  of  B ^ (0,  Qq): 

We  begin  with  the  asymptotic  representation  of  Legendre  functions: 

L£p!c .!*<«»]  (B-l) 


where 

41  = Vs  -t  ( i - 2m)i 
Equation  (13)  then  becomes 


( SIM0o 

±Jk,  cost1  _ 

** 

■s  + on*  ro+  0 SjnJjL 

Urn©  / 

00  [ COS  if. 

8 cos  4*p 

(B-2) 


First-order  and  second-order  roots  of  P mi  (cos  0) 
as  follows: 

First-order  Hoots 


0 arc  then  found 


sm  yf0 


to 

" rr 


= <"  -i  ♦?)£  • 


n 


1,2.3, 


(b-3) 


Second-order  Roots 


5m  if 


1 ~ ^’rri  co't®  COS  lf0 

8 v 


I Qn 


% - 


I - » 

8 «!" 


CO 


te„ 


- t 


By  taking  a Taylor  expansion  of  tan  about  (//q 
we  find 


(1) 


= nil- , 


, t*)  lo) 

4.  = ^ 


e. 


The  desired  expansions  of  and  — — which  appear  in  B'  (0,  9 ) 

1 cos  i//  cos  j)  in  o 


(B-4) 


(B-5) 


in  terms  of  first-order  and  second-order  roots  are  found  by  taking 

ip  r ip  _ dr  H ip  . Then. 

^ * Yo  To  > 
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iff  - % = i>B  +(i-Zm)1L  - tx80 

= rrie-o*)  + 


(B-6) 


H*o  = * « 1 6 

COS  i|»  = cos(if -H'o)  C0S4<0  - SMq(H'-%) 

- (-I>n  [cosi/'V-e.)  - e|e  5ii^] 

where  we  have  taken  sinejsc  and  cost  ss  ^since 
Similarly,  we  find 

sin  4;  - sin(iy-4>o)  co$4>0  cos^-iy,,)  sin^o 

= (-<jn  [ sin  l/°\e-9»)  ■¥  € - COS  i<,'0,(e-6.)] 


(B-7) 


(B-h) 


cos 4*0  = cos[  n n -r  e ] - (_|J 

With  (B-8),  (B-9),  and  ( B- 1 0)  in  (B-2)  and  the  asymptotic  representa- 
tion of  the  Bessel  function  the  terms  in  the  residue  series  approach  the 


(B-9) 
(B-  10) 


following  values  in  terms  of  the  first-order  and  second-order  roots, 
\ Thus, 

± — tOe,0e)  fUn  i(*o)  (zrrxaf  e.  %) 


V (o)  and  V (1  \ Thus, 
n n 


- i 


u, 

u 


± [ *■"  g*-)  i J±  I cos^W(e-e.j 

M i Sm9/  8,  ' 


-L  II  S'«  ^ (»-®o 
^ro>  1 1 


J 1 


-2v„ 


n » ) 


(B-ll) 
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1 

I 

s 

f 

I 


u 

I! 

SI 

If 


But, 


-*4. 


7 


whe  re 


-»o 


(o) 


e 


i 


] 


-fP, 1 

e 


[ . - Ok  7 

*r» 


ri  = IZJtm1  JL  col  9 . 

'Z  * 


Substituting  (B-12)  in  ( B — 1 1 ) and  retaining  the  zeroth-order  and  first- 
order  terms  in  1 /i'n»  the  right-hand  side  of  ( B - 1 1 ) becomes 


I 

9o 


/ sin  9P  V2  [ cos //“Ve-et)  e 

\ iin  0 / / 


Af  n -7*n  Cos^%,-ej7 

nn  L 1 12 


. ,0>; 


17  »/ 


The  first  term  in  the  braces  gives  0^(9,  9Q)  while  the  last  two 

terms  give  R (0,9  ) after  replacing  the  coefficient  1 by  ^o 
" ° i>  ( o ) n TT 


n 


Expansion  of  0^(0,  9q) 

On  using  (B-l)  in  (13)  of  Section  2,  and  dropping  higher  order  terms 
in  1 / V , we  find 


w 


sin  0Q  \ \ sirup 


sin  9 J 9q  cos  i p0 


?.r  . tan  i pn 

(l-4m  )[cot9  cot  i p -I  cot9  tan  ib  - 4(1 + 2m)  D 
• ( 1 _ O o 1 o'  ' wo 

8 v 


(B-12) 


(B-  13) 


-141 
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Substituting  (B-9)  and  (B-10)  for  sirup  and  cos  ip°,  respectively, 


B,  (0.  9 ) ^ 
]n  o 


/ sin  0Q\ 

1 

2 1 

sin  in  ^ ° ^ ( 9 - 0o)  + (4m2-l)  cos  v n ( 9 -Oo) 

lsin9  j 

9 

o 

rr ^ — 

l n J 

(B  - 1 5) 


If  we  now  use  the  asymptotic  representation  of  the  Bessel  function,  we 
can  express  the  nth  term  of  the  residue  series  in  terms  of  the  first- 


order  roots,  v resulting  in 


"in*  <0’S„ 


EL  B.  (9,  9 )i 
2 - ..  In'  o' 


i'  -1/4 
n 


vr\  . A -i  (ka-  tt/4) 

J"  n ka)  <2  e 


li 


('sin9o\ 
\sin9  / 


Isin  v (o)  (9-9  ) e_i 
9- 


(B  - 16 ) 


(o) 

n « 


where  v _i?iL  ^as  Been  used  in  the  magnitude  factor.  The  above 

8 


expression  may  be  summed  to  give  an  approximate  closed  form,  Rj^(9,  0^), 
as  follows: 


Ru(“'eo>-ERm^«o»  -Z{^nr) 

n = l nsl  ' ' 


sin  v ^ (6-9  ) e 1 1 n 77 
n o 

n tt 


2 tt  i 


/sin  9 X"1  °°  , f -i  v tt  - 0 

( O ) He  n O 

\ sin  0 / ni 


+ e)  -i"n(o)(*+  «o-0)- 

' - e 


with  gj(  a , ( ) , g^(  (1  , £ ) defined  below  equation  (47-c). 


( B - 1 7 ) 


i Hi 
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J 


Expansion  of  ( 0 , 0 Q) 

On  using  (B- 1)  in  (13a)  of  Section  (2)  and  dropping  higher  order  terms 


in  1 / g , we  find 


1S2n<  ° 


, l sin9o 

’ °o*  \ sin  » / 


COS 


0 


0 sin  ii/ 

O r o 


cot  (i/ 


o 


1- 


(344m2)  cot  0tan0  4 cot  0 Qcot  0 q 4 8(m-v)  0, 


8 g 


n 


We  now  find  the  first-order  and  second-order  roots  of  ( d / ri  0 q)  P ™ _ , 
= 0,  by  differentiating  (B-l)  and  setting  the  result  equal  to  zero.  Thus, 


u cos  dj 

r V-  O 


l-4m2cos  ip ^ 

2 "~ 
sin  0 8 )i 

o 


344m2 

8 


cot  0 sin  p 
o o 


0 


o r 

cot  0 = 

o 


l-4m2 

8g2sin20  . 

o 


I 


f 


First-order  root: 

Setting  < = 0,  cot  0Q  = O when  0 q(1)  = ( 2 n 41)  | , giving 

g = (n  4 J_  + m ) JT_  n = 1,  2,  ... 

n 4 2 0 

o 


(B-18) 


(cos  0q) 


(B-19) 


(B - 20) 


(B-  21) 
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Second-order  root: 


(1) 


By  expanding  cot  \p  in  a Taylor  series  about  t p ^ , we  find  from  (B-19) 


Cot0(1»  - cscV1'  (^(2)  - *0(1)>  =* 

o o o o 


(B-  22) 


and  (//  ^ = 0 ^ + f' 

o o 


(1) 


since  cot  ip  w— ► 0.  Then, 
o 


n 0 ^ + (1- 2m)  = (2n  + 1)  "2  + * 


n o 


(2)  (0)  . f 

or  M n = M n + ~ 


(B  - 23) 


To  find  expansions  of  — r 


T2-s-_.0-  ancj  sin_J^_  in  terms  Qf  first-order 


sin  0 sin  ii 

*o  ro 


and  second-order  roots,  again  take  p - ip  - p q + p Q.  Then, 


0 -0O  = Me  - eo) 


(%(o)+  it)  (e-  Bo> 


= M(o)(e-e)+-^-  ( 9 - 9 ) 

n o 9 o 

o 


(B  - 24) 


0 


( 2n  + 1)  -y-  + e = k — — + t , k=l,3,  ...  (B-25) 


Using  (B-  24),  (B-25)  and  the  fact  that  e and  e'  are  small,  we  find 


cos  0 = 


• i 71  j • ( 

-sm kr-  < sm  m 
2 j n 


(o)  ' « - - ' J-  ' — cos  M (o)(  e - e )>  (B-  26) 


( 9 - 9 Q)  + c 


n 


sin  0 = sin  k — |cosgj°^(9-0o)  - e “ sin  g ^ ^ ( 9 - 9o)j 


sin  0 = sin  ((2n+l)  -y  + € ) « (-1) 

o £ 


n 


(B-  28) 
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Using  the  results  of  (B-26),  (B-27),  and  B-28,  we  find 


4 


n 


n 4 


B 2„<  6 - 8 oU 


- V- 


n 


_d_ 

dx 


<-E»*  J-n<x> 


(2T:ka)2e"^ka+?  ) 


(B-  29) 


/sin  eo  y 

\sin  6 J 


^ o)e 

r n o 


si n Mn(0)(  0 - dQ)  + 


(o) 


IE 7Vos%  10 -V 


n 


n 


Ti  _ (3  + 4m  ) cotd 
'3  ~ 8 


n » 1 


We  also  take 
-2  H 


n 


where 


•i('I1So)+  t)7T 

e o 


-i  u (°K  i 

-i 

" [“7^] 


(B  - 30) 


V, 


8 6 


(3  + 4m  ) cot  6 . 


Substituting  (B-30)  in  (B-29)  and  retaining  only  the  two  lowest-order 


terms  in  — — ^ , we  have  for  the  right-hand  side  of  (B-27) 


1 


(o) 

-ifl 


n ff/Bingo\2  JL 
\sin  6 j eQ 


sin  fi^o)(e-dQ) 


(B  - 31) 


1— 2 p?3  coSfi^o)(0-eo)  - X v4  sin^o)(e-do) 

(oh  L 


( fj.  (0)) 

n 


n»l 


After  replacing  the  coefficient  1/^  ° by  0 Q/mr  , the  first  term 

in  the  braces  is  sumecl  to  give  a closed  form  approximation, 

R.  , (6,6  ),  while  the  last  two  terms  are  neglected.  Thus, 

Z </>  o 


QO  00  . V 

<i8’eo)='ZlRZn^e’(,J 

n= 1 n= 1 ' ' 

)i  oo  r ( n ' 

Vl  '<”■ 

/ jn  |_e 


(o) 

1 / \ -X/i  TX 

2 sin/u^°T0-  0Q)e 


^°'(7T-0  + 0)  -i/x  (°V+0  -6) 
' o n o 

- e 


i 

/ sin  0 \ 2 r 


(B-  32) 


where  g^a  f ),  g^/3,  £)  are  obtained  from  g^a,  ^),  g z(/3>0  by 
replacing  £ with  (£  + 1). 
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